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Abstract 


A 1-dimensional model for coherent quantum energy transfer through a complex of compressible 
boxes is investigated by numerical integration of the time-dependent Schrodinger equation. Energy 
is communicated from one box to the next by the resonant fluctuating Fermi pressure of the 
electrons in each box pushing on the walls and doing work on adjacent boxes. Parameters are 
chosen similar to the chain molecules of typical light harvesting complexes. For some parameter 
choices the system is found to have an instability leading to self-induced coherent energy transfer 
transparency. 
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I. INTRODUCTION 


Light harvesting complexes in plants and bacteria, as well as visual systems in animals, 
are typically arrays of chromophores organized spatially and held in place by proteins (see 
recent review articles (lJ, ^2 ] and references there). They absorb photons and transfer the 
excitation energy over distances as great as several nanometers with high quantum efficiency, 
estimated at 67% in the mammalian visual system, and up to 90% in the photosynthesis 
of higher plants. No artificial solar cell approaches these values. Biological systems achieve 
these values not just through being nanoscale, but also through being ultrafast, much faster 
than corresponding radiative lifetimes and other dissipative mechanisms. The isomerization 
of retinal in rhodopsin, the first step in vision, is one of the fastest photochemical reactions 
known, complete in 200 fs. A similar sub-picosecond timescale characterizes the first en¬ 
ergy transfers in photosynthesis. Understanding these fast energy pathways is a subject of 
intensive current research. 

In principle the chain of events following photoabsorption could be just incoherent tran¬ 
sitions from one excited state to another, passing energy to a reaction center (RC) in pho¬ 
tosynthesis, or to a new conformation of rhodopsin in vision. The speed of these processes, 
however, together with the experimental observation of coherences among some intermedi¬ 
ate states, has prompted the investigation of coherent, quantum mechanical evolution of the 
excited states as a contributing mechanism. These investigations have even raised funda¬ 
mental issues of quantum mechanics, suggesting that coherent energy transfer could be a 
kind of quantum computer, finding the best path by, in effect, trying all of them simulta¬ 
neously 3], or questioning whether coherent mechanisms seen in the laboratory have any 
relevance to performance under incoherent illumination by sunlight 4], The decay of coher¬ 
ence through interaction with the environment is relevant to such questions, and has been 
modeled by Hamiltonians that include interaction of the electronic state with ‘phonons,’ 
vibrational modes of the molecules, reviewed in 

This paper calls attention to an energy transfer mechanism that has not been part of the 
discussion. The excitonic state of a chromophore itself has a kind of phonon-like quality, 
even in the absence of nuclear motion, manifesting itself as an oscillating Fermi pressure. 
By steric interaction with a neighboring molecule, this oscillating pressure can do work and 
transfer energy. An essential feature of the model is that each molecule sees each neighboring 
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molecule as a wall, in the sense of classical mechanics, into which it cannot penetrate. (The 
reason for this is the Pauli principle, but the quantum mechanical origin of this excluded 
space is otherwise irrelevant to the molecular quantum states.) The motion of the wall, being 
a collective motion of the molecular electrons, will be essentially classical. In the model that 
we compute, the walls, considered as classical oscillators, have a natural frequency much 
lower than the exciton resonance frequency, another reason that their motion should be 
considered classical. The oscillating pressure, on the other hand, is entirely of quantum 
mechanical origin, as the beating of the molecular excited state with the ground state in a 
coherent superposition. 

Biological systems do not suggest any particular geometry for light harvesting antennae 
beyond the close proximity of chromophores and proteins. Accordingly we choose the sim¬ 
plest possible geometry to test the concept, in what is surely the simplest possible model of 
electronic excitation coupled to steric interaction in this way. It is exactly solvable in the 
time domain by numerical integration of the Schrodinger equation. The model is not pro¬ 
posed as a realistic geometry for an actual process, but only as a schematic way of thinking 
about a process that undoubtedly does occur at some level. We are not doing molecular 
modeling, but rather a much more primitive thing, investigating whether the proposed mech¬ 
anism could be significant under any circumstances, and might therefore be worth detailed 
consideration in more realistic geometries. 

The model is a chain of 1-dimensional boxes, each box representing a single chromophore. 
The relevant electronic states are quantum mechanical box states, where we allow the box 
to fluctuate in size, subject to the Fermi pressure of the electrons. The walls of the boxes 
represent the boundaries between chromophores, with masses assumed to be tens or hundreds 
of electron masses. In this way one box can push on the next, compressing it and raising its 
energy levels, or perhaps inducing quantum transitions. The mechanism is non-Coulombic, 
and universal in the sense that it is not constrained by selection rules. In particular it would 
work the same way for both optically allowed and optically forbidden states, if we included 
enough molecular detail to make that distinction. 

A computation shows that energy can be tranferred from box to box by this process, 
and in an unexpected way. As we show by computation, there is an instability in this 
system that allows the growth of fluctuations in Fermi pressure to feed on itself and grow 
to dominate the time development, creating a self-induced energy transfer transparency 
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through the system. For comparison we include also energy transfer by the familiar Forster 
resonance mechanism, and find that, within this model, the Fermi pressure mechanism can 
be of comparable importance. 


II. PARTICLE IN A COMPRESSIBLE BOX 

Consider a quantum mechanical particle of mass n in one dimension in a box that occupies 
the time-dependent interval [L(t), i?(i)] on the x-axis. This very elementary system has been 
treated before jb], but we describe it again here for completeness. The wave function $(a;, t) 
satisfies the time-dependent Schrodinger equation 


d$ h 2 d 2 $ 

lk ~dt ~ 

in the box [L, f?], as well as the boundary conditions 


( 1 ) 


$(L{t),t) = $(R(t),t)= 0. 


( 2 ) 


In terms of the variable 



(3) 


we can seek a solution to Eq. (JT]) in the form 


*(x,t) = *(£,<) 


(4) 


where T satishes 



1 U 2 <9 2 T 
(f?-L) 2 2^ d£^ 


(5) 


in the interval 0 < £ < 1, and the boundary conditions 


T(0,t) = T(l,t) = 0 


( 6 ) 


The wavefunction <f> determines the probability density |<F| 2 . In terms of the new variable 
£ the corresponding probability density is |T| 2 (f? — L). We verify unitarity in the form 



(7) 



( 8 ) 
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for all t, using Eq. (J5]) and several integrations by parts. 

For each t the wavefunction belongs to the Hilbert space of square integrable functions 
on the interval [0,1], vanishing at the endpoints. The natural inner product < | > on the 
Hilbert space at t is 

< f\g >— {R — l) /V^m (9) 

Jo 


and the functions 




y/2 sin mn£, 


are an orthonormal basis at each t. 

Expand the wavefunction T in this basis 


m — 1,2, 3,... 


( 10 ) 


'i' = y 

m 


(ii) 


Inserting this representation into Eq. (jS]) and resolving the terms in the basis {T m } leads to 
a representation of the Schrodinger equation as a system of ordinary differential equations 
for the amplitudes a m , 


d&m 

dt 


mn 


R — L + m 2 — n 2 

rijtm 


[L - R(- 1) 


n+ml 


ihm 2 7T 2 

(fi - Lj% a " 


( 12 ) 


Equivalently one can write this system in the “interaction representation,” defining 


Eq. (TT2]l becomes 


bm 

0(f) 


hi : 2 r* l 

- / 7 dt 

2 n Jo {R - L) 2 


db m 

dt 


2 


R-L 


y gi(m 2 —n 2 )(j> Tml H _ R(—l'\ n+m ] b , 

j m 2 — n 2 


In terms of these amplitudes the statement of unitarity takes the form 


(13) 

(14) 


(15) 


d * 

~TL Z_^ a m a m 

H 6 YY) 


0, 


(16) 


easily verified using Eq. (fT2]h and similarly for the analogous statement in terms of b m - 
The expectation value of the energy 


<E>= 


^ * h 2 m 2 n 2 
m a m a m 2^R — L) 2 


(17) 
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is time dependent because the moving walls can do work on the particle. In fact, using 
Eq. (H2D, we End 


d < E > 
dt 


2 9 


r 7 t 


jx(R — L) 3 

h \ 2 


J2mn[L - R(-l) n+m }M(a* n a m ) 


m,n 


E 


na r 


R 


£(-!)’ 


na r 


(18) 


(19) 


fi(R — L) 3 

The terms with n/m are rapidly fluctuating, while the terms with n = m simply reflect 
the adiabatic compression or expansion of the box. The coefficients of L and —R can 
be interpreted as the pressures at the left and right walls respectively. It is exactly the 
fluctuating component of these pressures that will be responsible for the energy transfer in 
the model that we compute. 


III. 2N ELECTRONS IN A BOX WITH MOVING WALLS 


The time development of multielectron states can be described by regarding Eqs. era 
and (fT5l) as equations for the evolution of fermionic operators (we develop this point in more 
detail in section |V|) . The ground state of a box with 2 N electrons would then be 

\N) = Ojv| a E'“ a it a i4-l^) (20) 

and one of the (doubly degenerate) first excited states would be 

\N ) = ajv + it«E- at i t «l|0), (21) 


where |0) is the empty box, and where the electron spin ('l' and J,) plays no essential role in 
what follows and will be largely ignored. Restricting attention to just the transition between 
these two states (HOMO-LUMO transition), and suppressing the subscript for spin, we have 
the time-dependent Schrodinger equation in the form 


dciTV —2 N(N + 1) 

dt R — L 2 N + 1 
dctTV+i 2 N(N + 1) 

dt R — L 2 N + 1 


(■L + R)un+i 
(L + R)cin — 


ihN 2 7T 2 


(R — L) 2 2/i 
2_2 


ClN 


ih(N + 1) 2 7T 

(R-L) 2 


-<Rv+i 


( 22 ) 

(23) 


or equivalently, in the interaction picture, 

db N _ -2 N(N+ 1) 

dt R — L 2 N + 1 
db]\f .)_i 2 N (N + 1) 

dt R — L 2 N + 1 


e -t(2jv+i)*^ + fy bN+1 
e i(2N+l >4>(L + R)b N 


(24) 

(25) 
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The time evolution of 0 is still given by Eq. (UTl) . If the box is squeezed symmetrically 
(L = — R ), no transition is excited, but the energy goes up. It is amusing to note that if the 
box is shaken rigidly (L = R) at the frequency 


Wres ~ (2IV + 1)0 


(21V + l)l«r 2 
2 ii(R-L) 2 


(26) 


Rabi oscillations are excited, i.e., oscillations in the amplitudes of the two-level system at a 
frequency proportional to the amplitude of the resonant driving force. 

The expectation value of the energy of the system is 


(E) 


n 2 

2 n(R - L) 2 


N(2N l + + VW| 2 + (N + 1)W,| 2 


and its rate of change due to the work done by the walls is 


d(E) 

dt 


P l L - PrR 


where the “pressures” are 


n V 

N(2N 2 + 1) 2 

L n(R — L) 3 

P R ^ 

2 i \\d\ \ ija_/v+i i ivajv| 

N(2N 2 + 1) 2 

R fi{R — L) 3 

2 1 | \ -LJO/v+l iVtt7v| 


(27) 


(28) 


(29) 

(30) 


Note that the fluctuating component of pressure is proportional to |ajvajv+i|, so that it 
is non-zero only in a superposition of excited and ground states. 


IV. 2N ELECTRONS IN EACH OF J BOXES 

If there are J boxes, the Hilbert space of states is the J-fold tensor product of the Hilbert 
space for a single box. The dynamics is just that of the formalism already described, but 
with subscripts on all quantities labeling which box they belong to. We assume no tunneling 
between boxes (i.e., the wave functions vanish at L and R as before). The dynamics takes 
place independently in each box, and does not, for example, lead to entanglement of the 
states, but the systems may be coupled through the motions of the walls. In particular, the 
case we shall consider, the boxes may be concatenated together, so that 0 = Li < R\ — 
L 2 = x\ < i ?2 = L 3 = X2 < ... < Rj-i = Lj = xj -1 < Rj = xj = J. Thus the J boxes 
occupy the interval [0, J] on the x-axis, and their average length is 1. 
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We also enlarge the dynamical system to include the moveable walls at x\, X 2 , ..., £j-i, 
while keeping xq = 0 and xj = J fixed. We treat the moveable walls classically, giving them 
masses M 1; ilf 2 , Mj_i sufficiently large, and imagine that the relevant forces on them are 
just the pressures due to the delocalized electrons of the previous section, i.e. 


dxj 

dt 

dvj 

dt 


= Vj 


= (Pr, ~ Pl m )/M. 


(31) 

(32) 


There is a conserved energy in this system, 


j-i 


E = £ <£)} + 2 E M jVj . 


3 (33) 

j =1 " j = 1 

useful for checking correctness of numerical computations. 

In light harvesting complexes the resonant frequencies (the energy scale) of the observable 
optical transitions typically decrease in the direction that the excitation energy follows to the 
RC. It was once suggested that this amounted to a kind of “energy funnel,” but the funnel 
is not unidirectional, according to more recent ideas: the photosynthetic complex may be 
more like a reservoir in which the captured energy is distributed js]. We can build such a 
structure into the chain of J boxes by choosing effective electron masses //i < /i 2 < ... < Hj 
that increase slightly in the direction that we expect the energy to flow. This alters the 
energy levels so that they are not initially in resonance, creating the situation that we had 
set out to investigate. 

Imagine that all boxes are in their quantum mechanical ground state. Mechanical equi¬ 
librium in Eqs. fl3T]) - ([32]) requires 

1/3 

j = 2,.., J. (34) 


x .i x .i-i 


dj -1 


x j —1 Xj —2 


dj 


Thus if the /ij’s increase with j, the boxes also become more compressed with j, in order to 
balance the pressure at each wall. Their transitions are not in resonance, however, because 
resonance between box j and box j — 1 requires 

1/2 

I , (35) 


Xj 


Xj -1 


dj -1 


dj 


x j —1 X j —2 

with box j slightly more compressed than it is in mechanical equilibrium. 

The following scenario motivated the computation. Taking the equilibrium state as the 
starting point, imagine that box 1 is suddenly placed into its excited state by absorption 







of a photon. The pressure is now higher in box 1, so that box 2 begins to be compressed. 
At sufficient compression box 2 comes into resonance with box 1, and its excited state is 
populated by the oscillating pressure through a Rabi transition. As the pressure in box 
2 increases, box 3 begins to be compressed, etc., passing the excitation along the chain. 
Computation shows that something like this happens, but also that the classical intuition is 
not completely correct. Rather, low frequency oscillations of the walls, induced by the initial 
photoabsorption, bring the boxes in and out of resonance. In the assumed pure starting 
states the pressure in each box is nearly constant, and transitions are driven only weakly, even 
during the resonances. As the states gradually become coherent superpositions of excited and 
ground states, an instability is reached in which the rapidly oscillating pressure, growing in 
amplitude, takes over the time development, and energy is then rapidly transferred through 
the whole complex. 


V. OPERATOR FORMALISM AND THE FORSTER TERM 


This section elaborates on the operator formalism suggested in section UTT1 The Fermionic 
operators a m and their adjoints obey canonical anticommutation relations 


Li $mm mi {®mi 0) 


(36) 


the empty box state obeys 

a m \ 0) = 0, (37) 

(and similarly for the operators b m of the interaction picture) for all m. A second subscript 
will indicate in which box the operator operates, as in a' m j, the creation operator for the 
mth state in the jth box. Operators for different boxes commute, since they operate on 
different factors in the tensor product. 

Hamiltonian operators can be written in these terms. The kinetic energy operator is 


J N +1 

Ho — E/ Hm,j a rn ,j a m,j 

spins j =1 m= 1 

Ti 2 n 2r m 2 

md = ^{R-L 3 r 


(38) 


(39) 


The effect of the moving walls on the time development can be given a Hamiltonian form 

J 2N(N T 1) 

Hi = ih E _ L .)( 2 N + 1) + E) { a N+U a N,j ~ alfjaN+ ij) (40) 
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a j — 


J 

ihJ2 ( a j b N+i,j b N,j - a*b]yjb N+ 1 j) , where 
-2- N(N + 1) ^ ^) iQN+ito 

Rj-Lj 21V+ 1 1 3 J) 


(41) 

(42) 


Finally we can include the Forster Hamiltonian, which couples nearest neighbor dipole 
moments, 


H, 


dipole—dipole 


Lj ' hj+1 
l r i - r j+il 3 


(43) 


The relevant part of the dipole moment of the jth box is 

e(N*\5t\N)j = e(R j —L j )2 £ sin(lV7r£) sin((N+l)n£) d£ = -e(Rj - Lj ) ^ . (44) 

Thus the Forster Hamiltonian is 
j -i 

Ih = ih £ [# i^N+ij+i^Nj+i ® &jvj6iv+ij A,- fejvj+i^w+ij+i ® ^iv+iy^Wj] j where (45) 


1=1 


Pi = e2 7i 


8N(N + 1) 
7i 2 (2N + l) 2 


Lj_|_i) (i?j Lj) 


(Rj +1 Ai 


3 i(2JV+l)(^ i+ i-0 i ) 


(46) 


and 7 j is a real phenomenological factor of order 1 which could have either sign, depending 
on the mutual orientation of the dipoles. 


VI. CHOICES OF PARAMETERS 


We follow Ref. [7] in modeling a typical light harvesting molecule as a box built out of 
2 N = 28 units (i.e., N = 14), each unit contributing 1 electron to the delocalized states, 
and each of length a = 0.1 nm. The resonant frequency is then 


U'k 2 {2N + 1) 

“ 2fj,[(2N — l)a] 2 


2.3 x 10 15 s" 1 


corresponding to a wavelength A = 830 nm. Choose units in which U = m e 
The unit of time is then 


(47) 

(2V-l)a = 1. 


(2 N - 1 ) 2 a 2 m e 
h 


63 fs 


if N = 14 


(48) 


We have made the conversion to physical time in reporting the course of the excitation 
through the chain of boxes. 

Forster transfer is the incoherent transfer of excitation by the mechanism of H 2 , usually 
calculated in time-dependent perturbation theory by Fermi’s golden rule. Over the short 
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times that we are modeling we instead consider H 2 to contribute to the coherent time 
development. We can choose parameters so that the coherent transfer time due to H 2 alone 
is a typical Fdrster transfer time, say 5 ps. Let there be just 2 boxes, rigid, each of length 
1, so that they are in resonance. H 2 then drives oscillations at the frequency 


LO 2 = Oi 


8 N(N + 1) 


4e 2 


71 ^ 


(49) 


ir 2 (2N + l) 2 

and this is also the coupling constant in H 2 . If we want the corresponding half period 
tt/co 2 to be 5 ps (physically), or 5/0.063 ~ 79 in our dimensionless units, then the coupling 
constant must be 

tt/79 « 0.04. (50) 


7T 1 


A common sense check of this parameter value comes from the virial theorem. The total 
kinetic energy of the 2N delocalized electrons is N(N + 1)(2 N + 1)7t 2 /6, and therefore the 
total electrostatic potential energy should be 


A/(iV + 1){ ™ +1)yl ~ = (2N e)(2N e nuc ) , (51) 

If we use the coupling constant from Eq. (I50lh and take the effective charge of the screened 
nuclei in each unit to be e nuc = e, we find, still using N = 14, 

' 1 \ N(N + 1){2N + 1) 


26 


(52) 


3Af 2 ir 2 (0.04) 

as if the chain molecule had length 1, by choice of units, but width only 1/2N, a reasonable 
picture. 


VII. ALGORITHMS AND RESULTS 

We integrate the Schrodinger equation 

ih^ r = (H 1 + H 2 )^ (53) 

at 

in the interaction picture, described above by the operators b m , arguing that over the short 
times that we will investigate, H 2 should also contribute to the coherent time development 
of the quantum state of the system. The Hamiltonian H 2 entangles the box states, so that 
it is no longer possible to treat the dynamics in each box separately. Where the state space 
could have dimension 2 J in section HVl it now must have dimension 2 J , a notable increase in 
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complexity. Let us choose a basis consisting of tensor products of J box states, each being 
one of | N) and | N*), with the Jth box represented at the left and the 1st box at the right, 
and let us label them by integers 0,1,..., 2 J — 1. The labels are read as follows: express the 
label in binary, and interpret the digits 1 or 0, left to right, as meaning the excited state or 
the ground state. Thus the label 13 = 1101 2 is for the state \N*) <8 \N*) 8 | N) 8 \N*). 

It is now straightforward to write out the matrix of H\ + H 2 in the interaction picture. 
For J = 2 it is 


h 1 + h 2 = in 


0 

—a* 

-a* 2 

0 

Ctl 

0 

-A 

-a* 2 

01.2 

A 

0 

—a\ 

0 

Oi 2 

Ctl 

0 


(54) 


More generally one can prove the following inductive scheme for constructing the matrix of 
H\ + H 2 (apart from the factor ih) for any J. Let Aj be the 2 J x 2 J matrix for (Hi + H 2 )/ih 
in the case of J boxes. Let Bj be the 2 J x 2 J matrix given by (in Matlab notation) 


Bj = diag(a;j +1 * ones(l, 2 J )) + diag(/3j * ones(l, 2 J 1 ),2 J x ) (55) 


i.e., aj + 1 on the diagonal, and j3j on a super diagonal. Then 


Aj + 1 — 


^4, 

B, 


-B) 

At 


(56) 


The pressures P/\ ; and Pr :i in box j are still given by Eqs. [2911301 using the variables of 
box j, where it is now understood that we must trace over all the other boxes j' ^ j. An 
efficient way to compute these pressures is to find the coefficients of Lj and —Rj in 


Thus, for example, if J 


d(E) UH 0 \ 1 

— -\-wj + wS [Ho ' Hi]) 

2 and the wavefunction is 


vT; — r T 

F v °mn ^ mn 5 

m,n =0 


(57) 


(58) 


using the binary notation for the labels, 

hV \n(2N 2 + i) 


P L1 = 
PR, = 


^(Ri-Lif 

h 2 7T 2 

m(Ri-Li ) 3 


+ \Nc 00 + (N + l)c 0 i| 2 + I Ncio + (N + l)c n | 2 


—-— -- + |A^Coo — (N + l)c 0 i | 2 + \Nc w — (N + l)cn| 2 


(59) 

(60) 
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We have investigated chains of 4 boxes, with slight systematic trends in their resonance 
frequencies from one box to the next, parameterized by systematic trends in the effective 
electron masses 

/ij = f + ;jAji j = 0,1,2,3 (61) 

where the detuning parameter A/i takes values 0.012 < A/i < 0.04. 

We first describe the time evolution without the Forster term (i.e., /3j = 0). A typical time 
sequence is shown in Fig. CD Here the detuning parameter A/i = 0.02, so that the effective 
electron masses were 1,1.02,1.04,1.06, and the wall masses were all M = 100. Contrary to 
intuition, the excitation placed initially in box 1 does not then move to box 2. Rather, after 
a delay, it suddenly appears in box 4. Thereafter, since energy is conserved, and there is 
no dissipation mechanism, it bounces back to box 3, back to box 2, and eventually becomes 
more chaotically distributed, but tending to stay at the bottom of the funnel (box 4). The 
reason for this is that once the states evolve into superposition states (and this takes some 
time, around 1 ps in Fig. [[]), there is no barrier to energy flowing coherently through the 
complex, as the oscillating Fermi pressure dominates the time evolution, inducing quantum 
transitions throughout the system. The change in the nature of the oscillations can be seen 
in the velocity of the third wall, Fig. [21 The slow oscillation of the wall about its equilibrium, 
governed mainly by the static steric pressure, gives way to rapid oscillations driven by the 
fluctuating component of Fermi pressure once the fluctuating component is large enough, 
at around t — 1.3 ps. The existence of a fluctuating component in one box drives Rabi 
transitions in the neighboring boxes that increase their fluctuating components, so that the 
process feeds on itself. This is the reason for the instability. It would be exponential growth 
if it were not bounded by unitarity. 

If we change the initial conditions so that box 1 is initially in the state x/0.991 N* > 
+v / 0l)I|jV >, i.e., we add a slight admixture of the box 1 ground state, the sequence 
that required a 1 ps delay now happens immediately, as in Fig. [3j This demonstrates the 
existence of the instability by starting the system a little further along in its unstable time 
development. In either case energy transfers, when they occur, are on a time scale of about 
a hundred femtoseconds. 

For each value of the detuning parameter A/i there is a range of wall masses M (all chosen 
the same, for simplicity) that allows this self-induced transparency. Despite the apparent 
universality of this mechanism, it still requires parameters to be chosen appropriately. If 
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M is too large, neighboring boxes are never brought into resonance by the oscillations that 
follow the initial photoabsorption. If M is too small, the classical oscillations are large and 
neighboring boxes spend too little time in resonance to evolve into the mixed states that 
drive the Fermi pressure mechanism. 

Now we add the Forster term to the coherent evolution. Using the coupling constant 
estimated in Section [V] we find the time evolution of Fig. [4l only slightly different from 
Fig. |H as one might expect, given the very different characteristic times for the two transfer 
processes that we are considering. With a stronger Forster coupling constant than we 
estimated above (0.1 instead of 0.04), the transfer is very noticeably affected, and the transfer 
is now mainly to the adjacent box 2 and not to box 4, as seen in Fig. 0 The two mechanisms 
are not exactly in competition. The coherent Forster interaction, by mixing excited states 
and ground states, actually speeds up the onset of the Fermi pressure process, as one sees 
in the peak corresponding to box 4, now just beyond 1 ps. 

VIII. DISCUSSION 

The Fermi pressure mechanism that we have described does transfer energy down a 
detuning gradient. The sudden transfer of energy from the first to the last box in the model 
is a counterintuitive and surprising feature. It is due to a self-induced transparency that is 
initiated by an instability in the oscillating Fermi pressure amplitude, leading to its rapid 
growth. This phenomenon could conceivably offer a dramatic speedup in energy transfer 
within close packed molecular complexes. 
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1.4 



FIG. 1: The excitation |a. tv+i | 2 as a function of time in box 1 (green), box 2 (red), box3 (cyan) and 
box 4 (blue). Here N=14, and the effective masses, creating the funnel, were /i\ = 1, H 2 = 1.02, 
/i 3 = 1.04, and 114 = 1.06. All wall masses were M=100. The excitation reaches the last (4th) box 
directly from the first box, with 50% strength, at about 1.3 ps. The Forster term is not included. 
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FIG. 2: The velocity of the wall between box 3 and box 4 around the time of the transparency insta¬ 
bility shows the slow, low amplitude oscillations giving way to the fast, large amplitude oscillations 
driven by the growing fluctuating component of the Fermi pressure. 
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FIG. 3: The excitation |a_/v+i| 2 as a function of time, as in Fig. [H but with a slight admixture of 
the ground state in box 1, so that the Fermi pressure fluctuations were already appreciable at t=0. 
The excitation now reaches the last (4th) box directly from the first box, with about 45% strength, 
in less than 400 fs. Transparency of the system is clear as the energy moves rapidly thereafter from 
box to box. 
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FIG. 4: The excitation |a. tv+i | 2 as a function of time in box 1 (green), box 2 (red), box3 (cyan) and 
box 4 (blue). Here N=14, and the effective masses, creating the funnel, were yiii = 1, /X 2 = 1.02, 
/i 3 = 1.04, and /x 4 = 1.06. All wall masses were M=100. The Forster term is included with a 
dimensionless coupling constant 0.04, as estimated in the text. 
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FIG. 5: The excitation |ajv+i | 2 as a function of time in box 1 (green), box 2 (red), box3 (cyan) 
and box 4 (blue). As before, N=14, and the effective masses, creating the funnel, were m = 1, 
^2 = 1-02, /i 3 = 1.04, and //4 = 1.06. All wall masses were M=100. The Forster term is included 
with a dimensionless coupling constant 0.1, corresponding to a Forster transition time comparable 
to the time for Fermi pressure oscillations to develop. The initial transfer is now mainly to box 2, 
not box 4. 
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